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We study the effect of freezing the positions of a fraction c of par- 
ticles from an equilibrium configuration of a supercooled liquid at a 
temperature T. We show that within the random first-order tran- 
sition theory pinning particles leads to an ideal glass transition for 
a critical fraction c = cr-(T) even for moderate super-cooling, e.g. 
close to the Mode-Coupling transition temperature. We first derive 
the phase diagram in the T — c plane by mean field approximations. 
Then, by applying a real-space renormalization group method, we 
obtain the critical properties for c — ck(T)\ — > 0, in particular the 
divergence of length and time scales. These are dominated by two 
zero-temperature fixed points. We also show that for c = cpc(T) 
the typical distance between frozen particles is related to the static 
point-to-set lengthscale of the unconstrained liquid. We discuss 
what are the main differences when particles are frozen in other 
geometries and not from an equilibrium configuration. Finally, we 
explain why the glass transition induced by freezing particles pro- 
vides a new and very promising avenue of research to probe the 
glassy state and ascertain, or disprove, the validity of the theories 
of the glass transition. 

Almost any liquid shows an impressive growth of the re- 
laxation time when supercooled below the melting point: 
indeed typical time-scales increase from picoseconds to hours 
in a rather restricted temperature window. This increase is 
so steep that below a certain temperature, T g , it is not possi- 
ble to equilibrate the system, which then freezes in an amor- 
phous solid, called glass pQ. The experimental glass transition 
temperature T g is not a critical temperature at all, but only 
defines a cross-over between the system relaxation time and 
the observation timescale. A recurrent question in the liter- 
ature is whether a true thermodynamic phase transition to 
an ideal glass state takes place at a temperature Tk < T g . 
This would explain the fast growth of the relaxation time and 
several other phenomena observed experimentally. Other sce- 
narii, based on purely dynamical transitions at zero (or finite) 
temperature or even no transition at all, might also provide a 
correct explanation [2j[3]. 

The study of the glass transition is hampered by the slowing 
down of the dynamics which is much more severe than for 
usual critical phenomena. Recent work on growing length- 
scales [4j suggests that even if a transition is present the cor- 
responding critical length that can be observed in experiments 
is no longer than a few inter-particle distances, meaning that 
experimental systems can be equilibrated only rather far from 
the putative transition point. This makes quite difficult to 
prove (or disprove) the existence of the transition and also to 
contrast theories solely on the basis of their critical properties. 
Here we propose a way to bypass these problems: we show 
that an ideal glass transitiorfj can be induced by a suitable 
random perturbation even for moderate supercooling. More- 
over, we show that the induced ideal glass phase can be ac- 
cessed and sampled easily. Contrary to the usual situation, 
where the existence (or the absence) of the ideal glass transi- 
tion can only be advocated on the basis of doubtful extrapo- 
lations, one can now approach it from both ends. As a conse- 
quence, determining whether there is indeed a true transition 
becomes feasible by using the usual machinery of standard 
phase transitions, in particular finite-size scaling. A crucial 
point is that this "high temperature" ideal glass transition is 
expected within the random first-order transition (RFOT) 5 
theory only. Thus, its study should lead to very sharp tests 
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allowing to ascertain the validity of the theory or disprove it. 
The suitable random perturbation we focus on consists in 
freezing the positions of an infinite subset of particles of an 
equilibrated configuration. This kind of procedure, first per- 
formed in [6], has recently been used to test the mosaic pic- 
ture of the RFOT theory 7 8, show the existence of medium 
range amorphous order [5] and study the growth of a static 
length [9] [10], called point-to-set [11]. In a very recent work 
Berthier and Kob have focused on several different geometries 
of frozen particles [12] , Their most remarkable finding is ob- 
tained when the positions of all particles outside two walls 
are frozen and when particles are pinned at random with a 
concentration c (see also [131 1 141 [15] ). In both these cases the 
relaxation time increases very fast when the fraction of frozen 
particles is increased or the distance between the two walls is 
diminished, even at temperatures close to the onset of glassy 
dynamics. The increase is so dramatic that it was claimed it 
could correspond to a true divergence [121 1131 115] . 
Here we show that within the RFOT approach this is indeed 
the case and that this divergence is related to an ideal glass 
transition. We first present phenomenological arguments sup- 
porting the presence of this phenomenon. We then theoret- 
ically derive the existence of this transition and its critical 
properties and discuss the growth of length and time scales by 
several methods. We make use of mean-field techniques devel- 
oped for disordered systems [16] and of the renormalization- 
group approach recently introduced in [17] . We show that the 
resulting theoretical problem appears to be related to nucle- 
ation close to a first-order phase transition in the presence 
of quenched disorder. The role of disorder is, actually, cru- 
cial and will be discussed thoroughly: disorder is expected to 
wipe out the transition in the case of two walls in dimensions 
three and less and, in the case of randomly frozen particles, 
in dimensions two and less. In the remaining cases, i.e. above 
the lower critical dimension, disorder should give rise to very 
unusual finite-size effects. Finally, we propose several ways to 
study these ideal glass transitions induced by freezing particles 
and we discuss the main reasons that make such studies fea- 
sible, in contrast with the usual case of the low-temperature 
ideal glass transition. We illustrate our main ideas by pre- 
senting numerical simulations of the Ising model when spins 
are pinned at random. 

General Arguments. As stressed previously, our results 
only hold within the RFOT theory, which explains the static 
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Ideal glass transition' usually refers to the transition, conjectured to take place at a finite tem- 
perature T K , towards an ideal glass state. We shall extend this name to all transitions induced 
by random pinning even though there are some important differences that we shall discuss in the 
following. 
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and dynamic properties of supercooled liquids in terms of the 
competition between the surface tension T and the configu- 
rational entropy density s c (T). The former is a measure of 
the extra free-energy cost paid when two different amorphous 
phases are in contact through a common surface. The lat- 
ter quantifies the multiplicity of amorphous phases in which 
the liquid can freeze. In mean-field models and approximate 
computations [IS] applied to finite dimensional liquids, s c is 
found to vanish at a finite temperature Tk- On a length- 
scale I (in spatial dimension d) the configurational entropy 
and surface tension respectively scale as TsJ d and Tl e (with 
< d — 1). In consequence, on lengthscales smaller than 
l PS = (T/Ts c ) 1/(d ~ e) it is preferable for the system to de- 
velop amorphous order, localizing in a given amorphous phase. 
Instead, for I > lps, it is better to break the amorphous or- 
der, pay the extra free-energy cost and gain entropy. This 
gives rise to the so-called mosaic state envisioned as an as- 
sembly of patches of length lps [5j[8]: a kind of micro-phase 
separated state where the number of phases is actually huge 

(proportional to e a " lps ). Recently, instantons computations 
confirmed the scenario depicted above, the existence of lps, 
defined as a point-to-set length, and its scaling with Ts c 

rjjrjsirji]. 

What are the consequences, within the RFOT theory, of freez- 
ing at random the positions of a fraction c of particles of an 
equilibrated configuration? Freezing actually means blocking 
the particle positions. The remaining particles then evolve 
in this frozen background. For this constrained liquid both 
the number of possible amorphous phases and the extra free- 
energy cost of the interfaces between them, T(T, c), change 
with c. The multiplicity of amorphous phases clearly is less 
than that of the unconstrained system. By pinning a particle 
one reduces the possible low free-energy configurations of the 
surrounding ones and, hence, diminishes the configurational 
entropy density by an amount Y (not to be confused with T). 
Thus, at least for small c, one finds s c (T, c) ~ s c (T) — cY(T). 
Generically, one expects s c (T,c) to be a monotonously de- 
creasing function of c, eventually vanishing for a critical 
fraction of frozen particles c K {T) ~ s c (T)/Y(T) < 1. Re- 
peating the previous argument about the competition be- 
tween configurational entropy and surface free energy cost, 
one then finds that the configurations of the constrained liq- 
uid are organized in a mosaic with larger tiles: the point- 
to-set length, l PS (T,c) = (T(T, c)/s c (T, c)) 1/(d " 9) , increases 
as a function of c — ck and diverges at the point ck(T) at 
which s c (T, c) becomes zero. Thus, we find that the condition 
s c (T,ck(T)) — 0, or equivalently s c (Tk(c),c) — 0, defines a 
line of ideal glass transitions. Its starting point, at c = 0, coin- 
cides with the ideal glass transition of the unperturbed liquid; 
its endpoint corresponds to the temperature at which T(T, c) 
vanishes. Within mean-field theory this happens at the onset 
temperature T 21, !22], see the top panel of Fig.l. By going 
beyond mean-field theory and by taking into account the ef- 
fect of the quenched disorder induced by the pinned particles, 
we find (see later) that the endpoint of the critical line actu- 
ally takes place at a temperature Th < T , see bottom panel 
of Fig.l. The value of Th depends on the ratio T^/Tr-, i.e. on 
the fragility of the liquid: when Td/TK decreases (more frag- 
ile liquids), Th/Td increases. Note that whether or not Tk (0) 
is different from zero, i.e. whether or not a finite tempera- 
ture ideal glass transition exists for the unperturbed liquid, 
is not crucial for the present discussion. In conclusion, we 
find that ideal glass transitions can be induced at tempera- 
tures between Tk and Th by increasing the fraction of frozen 
particles in the systenff]- The case corresponding to freezing 
all particles in the system except those between two infinite 



flat walls at distance d is similar and will be discussed at the 
end of the paper. We end this section with two important 
remarks on the newly found glass transitions. First, pinning 
particles can be considered as a very nonlinear way to probe 
the supercooled liquid state. In consequence, the existence of 
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Fig. 1. Top panel: Phase diagram for the p = 3-spin disordered model 
1221 . The continuous line corresponds to ideal glass transitions taking place at 
Tk( c )- The two dotted lines (from top to bottom) correspond respectively to 
the onset temperature T (c) and the mode coupling critical temperatures T^(c). 
Bottom Panel: Phase diagram obtained by the renormalization group approach for 
Td/TK — 1-4 (corresponding schematically to OTP). The variable on the x-axis 
is proportional to the concentration of pinned particles. The squares are the results 
obtained by RG for the ideal glass transition line Tk{c). The continuous line 
corresponds to the phenomenological result Ck cx: s c /Y. The endpoint of the 
transition line takes place at a temperature Th < Td in contrast to the mean field 
case. 

a glass transition for c — ck {T) reveals the presence of a new 
characteristic static lengthscale Ik(T) for the unconstrained 
liquid. This lengthscale corresponds to the critical mean 
distance between frozen particles: Ik(T) = \/c I i d . Since 
ck(T) oc s c (T), the new static lengthscale, Ik, is expected 
to diverge at T = Tk ■ It is interesting to remark that the di- 
vergence is milder than the one of lps'- Ik oc l/sV d oc l P d g S ^ d . 
The second very important issue is that, as we shall discuss 
later, the phase diagram of Fig.l is not expected within any 
theoretical approach other than the RFOT one. Thus, the 
new phenomena discussed in this work should allow one to 
test glass transition theories in a very stringent way. One 
could wonder, on the other hand, why showing the existence 
of ideal glass transition as a function of c is not as difficult as 
to show it as a function of temperature. In the usual case the 
main problem is that the low temperature phase is not known. 
In the frozen-particle case, this is not the case! Actually, the 
configuration C n f of the nonfrozen particles just after freezing 
out the others is automatically an equilibrium configuration 
for the new system [6] 1211 124| , even for c > ck- In this regime 
the nonfrozen particles just vibrate around the ideal glass con- 
figuration, which is now known: it corresponds to C n / and a 
lump of configurations around it. In consequence, one can 
now approach the transition from both ends (instead of mak- 
ing extrapolations from high temperature) and show whether 



Note that Y(T) is a microscopic configurational entropy loss and is not expected to vanish even 
at high temperature. 

^We are not interested in the regime corresponding to very large c where other physical effects, 
such as percolation of the frozen particles and Lorentz-gas physics, come into play I13II14I . 
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there is indeed a transition, study its critical properties, devise 
clever algorithm to equilibrate the system, etc. We postpone 
a more detailed discussion to the end of the paper. In the 
following we present theoretical arguments to corroborate the 
previous phenomenological ones. 




Fig. 2. Configurational entropy for the REM as a function of E for 
c = 0,c = C/^,c > Ck (from top to bottom). The equilibrium energy is 
stuck at the value _E e g(/3). For c < Ck the partition function is dominated 
by configurations with energy E eq ([3) (the straight line has a slope (3). For 
c > Ck the partition function is dominated by the contribution from the reference 
configuration; higher states lead to a sub-leading contribution. 

Phase Diagram and Mean Field Theory. An easy way 
to study within mean-field theory the effect of freezing de- 
grees of freedom is to focus on the Random Energy Model 
(REM) rjS] (see [H] for a similar computation). The REM 
is the simplest mean field model displaying an entropy van- 
ishing transition a la Kauzmann [T] and it is often used to 
understand some aspects of supercooled liquids. It is a spin 
model in which the energies of the configurations are inde- 
pendent and identically distributed (i.i.d.) Gaussian random 
variables of variance N, where N is the number of spins (such 
a scaling is needed to obtain a correct thermodynamic limit). 
The number of configurations with energies between E and 
E + dE reads Af(E) = 2^ exp[-~E 2 /{2N)]/V2^N. The equi- 
librium energy is given by the condition P — d In N(E) / dE, 
i.e. ft — —E eq (P)/N. The entropy S is simply equal to 
\nN(E eq ). The model has a transition at the temperature 
Pk = 1/Tk = V21n2 at which S{E eq ) vanishes. For T <T K 
the energy as a function of temperature does not vary any 
longer. Now, take a high temperature equilibrium configura- 
tion a and freeze a fraction c of its spins. We call this set Cf. 
In order to study the thermodynamics of the nonfrozen spins, 
one has to sum over all configurations having the spins in C / in 
the same state than in a. The energies of these configurations 
are also i.i.d. Gaussian variables. As a consequence, one can 
repeat the previous analysis except that now the number of 
possible configurations is reduced from 2 N to 2 JV ' 1_c - ) . Thus 
the energy distribution reads: 

M{E,c) = 2 JV(1 " C) e X p(-E 2 /(2N))/V2^N 

The relationship —E eq (P)/N = P remains valid but 
the total entropy is reduced by a factor —iVc In 2, i.e. 
S(E,c) = S(E,0) - ATcln2, see Fig.2. At the value c K = 
S(E eq , 0)/{N In 2) there is an entropy-vanishing transition. 
For c > ck the system is frozen in the configuration a since 
all the other ones have an energy that is higher and deter- 
mined by S(E, c) = 0, see Fig. 2. The analysis of this sim- 
ple model confirms our previous assumptions: the configura- 
tional entropy per spin is indeed reduced by a factor — Yc, 
here Y = /3 2 /[2(l — c)], and the low-temperature ideal glass is 
known: it is the high-temperature configuration a. 



An interesting additional information can be gained making 
the same analysis on models that display also a dynamical 
glass transition, akin to the Mode-Coupling Theory (MCT) 
transition. Following [27] we expect that since this transition 
is due to the emergence of metastable states, by freezing par- 
ticles one blocks unstable relaxation modes and transforms 
unstable states in stable ones, thus raising the MCT transi- 
tion temperature with respect to its c = value T^. Within 
mean field theory this is only possible up to the temperature, 
T D (c), where T(c, T) vanishes. We have analyzed the dynam- 
ics and the statics of the spherical p-spin model (details will 
be presented elsewhere [22], see also [23]) and found that in- 
deed Td and Tk are shifted upward by freezing particles. The 
lines Td(c) and Tk{c) join at the point at which both tran- 
sitions disappear, see the phase diagram for the p = 3 spin 
model in the top panel of Fig.l. This is in agreement with the 
results found previously by Semerjian and Ricci-Tersenghi in 
the context of combinatorial optimization problems [28] . 
Critical Properties and Renormalization Group Anal- 
ysis. In order to go beyond mean field theory and obtain the 
divergence of length and time scales one has to take into ac- 
count non-perturbative fluctuations. To this aim we use the 
approach recently developed in [17] and generalize it to the 
case of pinned particles. The analysis is quite cumbersome. 
In the following we shall only discuss the physical scenario 
that emerges from our computations; technical details will be 
presented elsewhere [22] . Let us first informally review a few 
important points related to the RFOT theory and the RG 
analysis of [17] ■ The order parameter of RFOT is the global 
overlap q(C r ,C) between a reference equilibrium configuration, 
C r , and another equilibrium one, C, that is constrained to be 
equal to C r on the boundary of the sample. The overlap q 
becomes different from zerqj at the transition because the 
boundary pinning field forces the configuration C to freeze in 
the same ideal glass state to which C r belongs. This transition 
is driven by the decrease of the configurational entropy, which 
plays the role of an ordering field favoring the zero overlap — 
liquid — phase. Including, in a renormalization group sense, 
non-perturbative fluctuations over the length-scale I one finds 
|17] that the renormalized liquid is characterized by a config- 
urational entropy s c (t), proportional to s c l d , and a term Y(l), 
akin to a surface tension in the replica field theory formula- 
tion, proportional to Yl d_1 . The latter is related to the sur- 
face configurational entropy loss between high and low overlap 
region^- The values of s c and Y are obtained from mean field 
theory. By integrating out non perturbative fluctuations un- 
til the lengthscale I, one finds that a system whose overlap is 
fixed to one on the boundary of a closed region of size I <g.lps 
thermodynamically prefers to remain in a high overlap state 
inside this regions, i.e. amorphous order can be stabilized 
and the number of possible amorphous states is s c l d below 
lps- The contrary happens for I 3> lps since Y(l) — > 0. No- 
tice that at the transition, when the configurational entropy 
vanishes, the point to set length diverges. 
From the RG point of view, the flow of the coupling con- 
stants is the one characteristic of zero-temperature disconti- 
nuity fixed points [29| . Indeed, one can make a useful analogy 
with the behavior of the ferromagnetic Ising model below T c in 
a negative magnetic field, see [3D]. The role of the Y and s c are 
played respectively by the coupling J and the negative mag- 



Depending on the normalization the overlap is zero or assumes a trivial value in the liquid phase. 
,J This allows us to clarify why we have chosen similar notations for T and Y. The reason is that 
they are related technically, they can be both associated to a surface free-energy cost (but the cor- 
responding free-energy is different) and physically, since the surface free energy mismatch between 
amorphous states, i.e. T, can be interpreted as the fundamental origin of the entropy loss due to 
pinning, i.e. Y. 
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netic field h: the negative magnetization (low overlap) state 
is favored by the field but the metastable one, characterized 
by positive magnetization (high overlap), can be stabilized on 
lengthscales smaller than lps by imposing favorable bound- 
ary conditions. Within this analogy the (T — Tk) — > limit 
corresponds to h — > 0. 

How does the scenario depicted above changes when one pins 
particles at random? By doing that, one forces the local over- 
lap to a high value in microscopic regions which are Poisson- 
distributed inside the sample. Close to all these regions, the 
system is subjected to an extra field of the order of Y, favoring 
the high-overlap state. This field opposes the global entropic 
driving force, which instead favors the low overlap state. In 
consequence, it is easy to grasp that by increasing the fraction 
of pinned particles one can induce a transition toward the high 
overlap state before the configurational entropy of the uncon- 
strained system vanishes, as found within mean-field theory. 
The RG analysis becomes more complicated than the one for 
c = because the local field induced by pinning particles is 
located in random positions. Because of this quenched disor- 
der, Y(l) and s c (l) become random variables and one has to 
follow the RG flow of their probability law. We have general- 
ized the analysis of [17] to this case and found that the flow 
is qualitatively identical to the one of the zero-temperature 
Random Field Ising model, as one would expect on the ba- 
sis of the analogy discussed above. Indeed, pinning particles 
translates into fixing a fraction c of randomly located spins in 
the up position and, hence, leads to an onsite bimodal distri- 
bution for the magnetic field, which is negative and equal to 
h with probability 1 — c and positive and larger than 2d J with 
probability c (a field larger than 2dJ fixes the spin in the up 
position on a hyper-cubic lattice in d dimensions). 
Let us first describe and discuss our results for temperatures 
slightly above Tk- In this case, by increasing c toward its 
critical value we find results very similar to the ones obtained 
approaching Tk at c = 0: the average value of the sur- 
face tension on the lengthscale I increases as Yl d ~ x , whereas 
the average value of the configurational entropy increases as 
k(c K (T) - c)l d oc s c (c,T)l d . The corresponding variances be- 
come much smaller than the averages, so that the random 
variables Y(l) and s c (l) almost do not fluctuate around their 
means. Y and k depend on T and c and are of the order of 
one, whereas s c (c,T) is the average configurational entropy 
density for a given value of c and T; it goes to zero linearly 
at ck{T). The length for which Y(l) becomes of the order of 
s a (l) is the point-to-set length, lps{c,T), for the pinned sys- 
tem (which is much larger than the one for the unconstrained 
liquid). For I 3> lps one finds that Y(l) — > and hence one 
cannot stabilize amorphous states. The scaling of lps with 
(ck(T) — c) is the usual one: 
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As already discussed in [17], the physical picture resulting 
from the RG analysis is that the renormalized system on the 
point-to-set lengthscale is like a liquid at its "onset temper- 
ature". The difference with a normal liquid are the typical 
value of the free energy contributions, which respectively are 
of the order of Yl P ~g and ~s c l PS , both very large compared 
to the temperature. This suggests that the typical relaxation 
time follows an Arrhenius law and, therefore, a generalized 
Adam-Gibbs relation 1311: 



Important differences with the c = 0, T = Tk case emerge 
when varying c at higher temperatures above Tk- Physically, 
the main consequence of fixing a higher reference temperature 
is to increase the value of the configurational entropy of the 
unconstrained liquid. Hence, the fraction of pinned particles 
needed to counterbalance this effect and to induce the tran- 
sition increases, see Fig.l. By doing so the variance of the 
random field, proportional to c (for small c), increases too. 
We find that above a certain temperature, Th, disorder fluc- 
tuations induced by pinned particles are so strong enough to 
make the glass transition disappear, in agreement with what 
was found for the RFIM when increasing the disorder [321133] . 
We studied the RG flow approaching Th- We found that the 
glass transition line ends in a second order transition point and 
that the long-distance physics is dominated by the RFIM zero- 
temperature critical fixed point. This could be expected [31] 
on the basis of the results [35] . We found that the means and 
the variances of Y(l) and s c (l) scale as l e , with 9 — 1.491(4) 
at T = Th and c = CK{Th) (to be compared to the one ob- 
tained for the RFIM by simulations 6 = 1.49(3)|52]). When 
T is close, but smaller than Th, the critical properties of the 
glass transition become quite complicated since the RG flow 
is first attracted by the RFIM zero temperature fixed point 
and approaches the standard first-order discontinuity fixed 
point for very large lengthscales only. Calling v the corre- 
lation length exponent for the RFIM (and using the stan- 
dard notation for the other exponents [321 133) ), we find that 
only for lps ^> l/(Tft —T) v the usual scaling laws discussed 
above are obtained. Instead, for lps >C l/(Th — T) v , i.e. for 



(T h - T) 



Ips oc 



v(d+^-2 v )/2 
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■C Cfe — c <C 0(1), one finds: 



log T OC 



{C K {T) _ c )29/(d-H7-2„) 



Because in d=3 the exponent 6 is very close to 1.5 and rj ~ 2?7, 
the relaxation time is expected to grow in this regime al- 
most according to the so called Vogel-Fulcher-Tamman law 
but with respect to the variable ck (T) — c instead of T — Tk , 
i.e logr oc (ck(T) — c) _1 . 

A final important outcome of the RG analysis concerns the 
T — c phase diagram. In particular we find that the phe- 
nomenological law ck(T) oc s c (T)/Y(T) is quite well obeyed 
for T < Th and that the temperature Th corresponding to the 
endpoint of the critical line is depressed with respect to its 
mean field value T . The phase diagram in the lower panefl 
of Fig.l corresponds to Td/TK — 1.4, a value typical of fragile 
liquids like OTP [19] (T d OTP ~ 285 A, Tg TP = 203if). In 
this case Th is very close to T^: Tj, ~ 280A'. Within our 
approximate treatment the decrease of Th with respect to T 
depends on the ratio of the mode coupling temperature to the 
Kauzmann temperature. Smaller values of Td/TK correspond 
to more fragile supercooled liquids and lead to an increase of 
Th/Td, that can reach values larger than one (and viceversa 
for more strong liquids). 

Other Ways of Freezing Particles and the Role of Dis- 
order. Pinning particles at random in space introduces a lot 
of quenched disorder and, hence, leads to a substantial de- 
crease of the endpoint of the ck(T) line. In order to reduce 
this effect, one can pin the particles closest to an ordered 
template, e.g. a periodic one [3H]. This leaves as source of 
quenched disorder, the fluctuations of the reference equilib- 
rium configuration only [351 137| . In this case we expect sim- 
ilar results, the only change being the precise form of ck(T) 



^ \ 6 We cannot draw the equivalent of the MCT line because this is defined within mean field theory 

1o£ T OC OC OC only. The value of (at c — 0) that we use and draw in the lower panel of Fig. 1 is obtained 

T Ts^ - ^ (Cfc — c)^~ 1 ' from mean-field theory. 
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and the value of Th, which is expected to be higher. Actu- 
ally, the relationship with the RFIM Another case considered 
in the literature [12] corresponds to freezing all particles out- 
side two walls at distance I. Following the Kac analysis of 
Franz |21| . one finds qualitatively the mean field phase dia- 
gram of the top panel of Fig. 1, where now I is the control 
parameter playing the role of c . The transition takes place 
at a value of I proportional to the point to set length of the 
unconstrained systerrQ. The effect of disorder is not taken 
into account in the Kac analysis. Since the positions of the 
particles pinned beyond the walls are random, the constraint 
on the coarse-grained overlap is not homogeneous close to the 
walls and lead to an effective quenched disorder. Applying our 
RG treatment to the two wall case and including the effect of 
disorder, one would find results similar to the previous ones 
but with d — ¥ d — 1. This has two important consequences. 
First, as for randomly frozen particles, the mean field phase 
diagram is modified: the endpoint of the critical glass tran- 
sition line takes place at a temperature less than the onset 
temperature computed by mean field or Kac-like techniques. 
Second, because the system inside the two walls is d — 1 di- 
mensional, the lower critical dimension is higher than the one 
expected in the case of randomly frozen particles. No first 
order transition takes place in two dimensions for the Ising 
model in presence of quenched randomness [38]. Since within 
our framework glass transitions are related to the first order 
phase transitions of the RFIM, we conclude that glass tran- 
sitions can be only induced in systems with dimensions equal 
to four or higher in the case of two walls. Instead, in the case 
of randomly pinned particles, also three dimensional systems 
display glass transitions when c is increased. 
The final case we shall discuss corresponds to particles that 
are frozen not from an equilibrium configuration but instead 
completely at random, i.e. from an infinite-temperature con- 
figuration. This is akin to studying liquids in presence of 
external quenched randomness [39] , a case relevant for porous 
media [24] . Because the quenched disorder is not tailored to 
an equilibrium configuration at temperature T, we do expect 
qualitative changes. In fact the phase diagram should be qual- 
itatively similar to those obtained in [391 124] . which are quite 
different from ours and resemble that of the p-spin spherical 
model in an external magnetic field: above a certain frac- 
tion c the glass transition becomes continuous (the Edwards- 
Anderson parameter has no discontinuous jump) but does not 
disappear; moreover, the glass phase for c > ck cannot be 
reached from high temperature without crossing a transition 
line, contrary to what is shown in Fig. 1. We shall discuss 
our phase diagram in more detail in the following. 
Glass Transitions Induced by Random Freezing: Nov- 
elties, Consequences and Numerical Verifications. The 
most striking feature of the phase diagram shown in Fig. 1 
is that it is possible to follow a path that connects the liquid 
and the glass phase without crossing any phase transition, as 
for the liquid-gas phase transition^ (see [22] for more details). 
Afterall, this is not so surprising since the ideal glass phase 
coincides up to small fluctuations with the initial configura- 
tion from which one freezes particles. Thus, one does not have 
to equilibrate the system on impossibly large time-scales or to 
solve a daunting optimization problem in order to reach the 
ideal glass. This should allow one to perform much more de- 
tailed investigations of the glass transition. In particular, the 
belief that it is impossible to prove (or disprove) the existence 
of the ideal glass transition because the time-scale for equili- 
bration diverges so rapidly maybe ill-founded. For instance, 
the RFIM also has a transition with activated dynamic scal- 
ing and despite numerical difficulties several properties of the 
transition have been obtained by numerical simulations. The 



main difference between the transition of unconstrained super- 
cooled liquids and the RFIM's one is that for the former the 
low temperature phase — the ideal glass — is not known. This 
prevents one from approaching the glass transition from both 
ends, from devising clever algorithms to equilibrate the sys- 
tems, etc. In the case of frozen particles, this is no more the 
case. Thus, one can endeavor to demonstrate the existence 
of the ideal glass transition by applying standard statistical 
physics techniques. Ideally, one could first study the transi- 
tion at a finite temperature well above Tk , at a finite ck(T). 
This would be already a major result in itself. Afterwards, by 
decreasing the temperature, one could measure the ck(T) line 
and show, at least by extrapolations, the existence of the ideal 
glass transition for usual, unconstrained, supercooled liquids. 
Even though not impossible, numerical verifications will be 
very challenging (work is underway [36]). We think that one 
should start first by studying finite-size scaling. If there is in- 
deed a transition, the averaged probability distribution func- 
tion of the overlap with the reference configuration should 
behave as in first order phase transitions: a single peak at 
small (large) overlap below (above) ck in the thermodynamic 
limit. Only at ck, a two peak structure should be present. 
The corresponding finite size scaling is expected to be simi- 
lar to the one of first order phase transitions in presence of 
quenched randomness, see e.g. [40]. The effect of finite size 
on the relaxation time r is instead different. Actually, already 
the thermodynamic limiting behavior is different: r, measured 
for instance by the intermediate scattering function, diverges 
approaching the critical line from the liquid side but not from 
the glass sidfl Finite sizes make r finite at the transition. 
Coming from the glass side r should increase very fast on a 
window that shrinks to zero in the thermodynamic limit. 
Actually, an instructive case to grasp the main ideas and diffi- 
culties is provided by the three dimensional Ising model below 
T c . Take an equilibrium configuration in presence of a posi- 
tive magnetic field h, freeze a fraction c of spins at random 
and revert the field to —h. Now the original configuration 
is metastable for c smaller than a critical value c*, and sta- 
ble for larger ones. Despite some important differences, the 
situation is very reminiscent of the one we discussed for su- 
percooled liquids. Indeed, our RG analysis is precisely related 
to this problem at zero temperature with h playing the role 
of the configurational entropy, the positive and negative mag- 
netization playing the role of high and low overlap. The life 
time of the metastable state, which diverges for c* — c — > 0, 
is analogous, mutatis mutandis, to the relaxation time that 
takes to the liquid to decorrelate from the initial configura- 
tion and that diverges for ck — c — > 0. We have studied the 
behavior of the overlap with the reference configuration as a 
function of time, see Fig. 3. We also explicitely verified that 
the overlap distribution behaves in the way described above, 
see the inset of Fig. 3. We found that by increasing the value 
of c the system spends more and more time in the metastable 
high-overlap state before equilibrating at low overlap. The 
timescale for relaxation at low overlap diverges at a critical 
value c* (T, h) with the same law discussed before in the glass 



Instead, for randomly pinned particles, the distance between frozen particles at the transition scales 
as the point to set length at a power less than one, see the phenomenological arguments introduced 
before and 1411 . 

8 At variance with the liquid-gas case, changing temperature in the phase diagrams of Fig. 1 does not 
simply mean changing a control parameter. Different temperatures lead to different configurations 
of frozen particles, hence changing temperature means studying different systems. 
9 This notwithstanding, the self and the collective relaxation times can be both very large, although 
finite on the glass side. 

10 The major important difference between this simple Ising model and real liquids is that the orig- 
inal reference configuration is not an equilibrium one. The main consequence is that the system is 
not automatically at equilibrium after having frozen out particles; thus, relaxation and equilibration 
time are not the same, contrary to the case of super-cooled liquids. 
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transition case [3 logr oc l/(c — c*) d ~ 1 . The curves in Fig. 3 
reproduce qualitatively well the ones obtained in supercooled 
liquids by Berthier and Kob [12] (besides the trivial differ- 
ence that asymptotic low overlap value is negative). This is 
an indication that the RFOT theory, and the related physical 
mechanism we have unveiled in this work, provide a viable 
explanation for the mysterious and generic dramatic slowing 
down of the dynamics, which has been observed in simulation 
when pinning particles of super-cooled liquids |12l 1131 1141 115j . 
In order to show that this slowing down is indeed related to 
a phase transition, one needs to study finite size scaling as 
discussed above. 




1 10 100 1000 10000 100000 1e+06 

t (MC STEP) 



Fig. 3. Overlap with the reference configuration as a function of time for the 3D 
Ising model at T = 4, h = 0.12 on a cubic lattice with N = 8 * 10 3 spins (J is 
set equal to one). Increasing the fraction of frozen spins it takes a much larger time 
to escape from the high overlap phase. Inset: overlap distribution slightly above and 
below c* (T = 4, h = 0.09, N = 10 3 ). 

Conclusions. In this work, motivated by the numerical re- 
sults [121 1131 1141 115] we have studied the effect of pinning 
particles of an equilibrated supercooled liquid configuration. 
We have focused on geometries such that the number of un- 
frozen particles is infinite. In these settings, pinning can lead 
to a true phase transition and not just a cross-over, like in 
the cavity case studied previously [H [9]. Indeed, we have 
shown that pinning leads to a decrease of the configurational 
entropy and, within the RFOT theory, to an ideal glass tran- 
sition when the confinement parameter c is increased, even for 
relatively high temperatures, e.g. close to the Mode- Coupling 
transition temperature. The relaxation timescal43 is infinite 
along the ck(T) line and finite everywhere else, in particular 
for c > ck(T). Although in the latter regime motion of single 
particles can be quite slow, as it is in a crystal, we predict 
that if properly equilibrated the system should show a non- 
monotonous behavior of the relaxation time as a function of 
c for T < T h . 

Pinning particles may be expected to generically slow down 



the dynamics whatever is the correct underlying theory 
(RFOT theory, facilitation theories [2], frustration limited do- 
mains [3], etc.). However, the nonmonotonous behavior of the 
relaxation time and the concomitant thermodynamic phase 
transition are expected within RFOT onljQ, thus opening the 
way to new and crucial tests of glass transition theories. For 
example, in the case of kinetically constrained models (KCM), 
we expect that one can induce dynamical glass transitions by 
increasing c above a critical value c^. However, since the ki- 
netic constraint can be violated in reality, the relaxation of 
super-cooled liquids would proceed via alternative and very 
slow relaxation mechanisms for c > Cd and would lead to a 
behavior very different from the one expected within RFOT. 
A detailed study of the effect of pinning particles in Kineti- 
cally Constrained Models would be worthwhile. 
We notice that our approach shares important similarities 
with two other approaches: the one of Franz and Parisi [42], 
consisting in biasing the thermodynamics by favoring configu- 
rations having a high overlap with an equilibrium one; and the 
one of Garrahan and Chandler [44], consisting in biasing the 
dynamics by favoring dynamical trajectories with low dynam- 
ical activity. The mean-field result analysis of the former leads 
to a phase diagram that is indeed analogous to ours (the bias 
strength plays the role of c). We expect this relationship to 
persist even beyond mean field. Some numerical simulations of 
three dimensional liquids provide promising results in favor of 
this approach 42, 43 . Theoretical analysis of KCMs and nu- 
merical simulation of glass-forming liquids 44, 45] show that 
also in Garrahan and Chandler approach one finds a phase 
diagram similar to the one in Fig.l. There are however im- 
portant fundamental differences in this case. The main one 
is that the phase transition taking place at cjf(T) and shown 
in Fig.l exists within the RFOT theory only, contrary to the 
dynamical transition studied in 44, 45 that takes place in 
facilitation-based models but also in those described by the 
RFOT theory (see [46]). 

We end this paper by stressing again that one of the main 
advantages of the glass transition induced by random freezing 
is that the ideal glass phase can be obtained easily. Hence, if a 
transition indeed takes place at ck , it can be very likely stud- 
ied thoroughly, in particular by finite-size scaling. Since we 
find that the critical properties for ck — c — > + and T > Tk 
are the same ones as those of the ideal glass transition when 
T — > Tk for c — 0, studying the effect of freezing particles ap- 
pears to be a new and promising way to understand the glass 
transition of supercooled liquids. These studies can be per- 
formed in supercooled liquids by numerical simulations and 
also on colloids by experiments using optical trapping. 
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